Approximate Quantum Cloning with Nuclear Magnetic Resonance 
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Here we describe a Nuclear Magnetic Resonance (NMR) experiment that uses a three qubit NMR 
device to implement the one to two approximate quantum cloning network of Buzek et al. 
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Quantum information processing j^] has been the sub- 
ject of much recent interest, not only because it offers 
new modes of computation and communication, but also 
because quantum information differs from classical in- 
formation in several fundamental ways. One important 
example is the fact that it is impossible to accurately 
clone (copy) an unknown quantum state ^ , and so quan- 
tum bits (qubits) cannot be duplicated. It is, however, 
possible to prepare an approximate copy , and several 
schemes for optimal approximate cloning have been de- 
veloped. Nuclear Magnetic Resonance (NMR) |0, |, ||] 
has already been used to demonstrate simple quantum 
information processing methods jM pi H, |l^, |ll|, ^ 

|16H7[ S H il and here we de- 

scribe an NMR experiment that uses a three qubit NMR 
device to implement the one to two approximate cloning 
network of Buzek et al. 

We started by slightly modifying the approximate 
cloning network [ p5| of Buzek et al. to take advantage 
of our specific hardware; our version of the network is 
shown in Fig. ^ This takes in one qubit in an arbitrary 
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FIG. 1: Our slightly modified version of the approximate 
quantum cloning network developed by Buzek et al. The 
initial "preparation" stage has been replaced by an alter- 
native network, which is simpler to implement with NMR 
techniques; the second "copy" stage is unchanged. Filled cir- 
cles connected by control lines indicate controlled tt phase 
shift gates |2^] , empty circles indicate single qubit Hadamard 
gates, while grey circles indicate other single qubit rota- 
tions. The two rotation angles in the preparation stage are 
6*1 = arcsin (1/^3) ~ 35° and 02 = n/12 = 15°. 



FIG. 2: The three qubit system provided by E-(2-chIoro- 
ethenyl)phosphonic acid dissolved in D2O and its NMR 
spectrum. The broad peak near —50 Hz is a folded signal 
arising from residual HOD. 



state (normally considered to be a pure state, E^nd 
two ancilla qubits in state |0); the two ancilla qubits are 
prepared into an appropriate initial state, and then the 
input qubit is copied. At the end of the cloning sequence 
the three qubits are all entangled with one another, and 
so the reduced density operator descriptions of each qubit 
(obtained by tracing out the other two qubits) correspond 
to mixed states. The two ancilla qubits are both in the 
state + (that is, approximate clones 

of Itjj) with fidelity 5/6), while the input qubit is now 
the approximate transpose of \ip){ip\. Alternatively, us- 
ing the identity IV'} (V" I + IV'''') (V""^ I — li the states of the 
ancilla qubits can be written as ||V')(V'I + ^1/2- Cloning 
is not, of course, confined to pure states, and any mixed 
state p can be cloned to give + |l/2 . Such mixed 
state cloning is particularly easy to study using NMR 
techniques, as the ensemble averaging inherent in NMR 
experiments is a help rather than a hindrance. 

We implemented this network using a 3 qubit NMR 
quantum computer based on the single "^^P nucleus (P) 
and the two nuclei {A and B) in E-(2-chloroethenyl)- 
phosphonic acid (Fig. ||) dissolved in D2O. Ethynylphos- 
phonic acid diethyl ester was prepared from trimethyl- 
silylacetylene by successive treatment with EtMgBr 
and diethyl chlorophosphate, followed by cold aqueous 
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TABLE I: Measured NMR parameters for E-(2-chloro- 
ethenyl)phosphonic acid dissolved in D2O (0.08M, 20°). Fre- 
quencies are measured with respect to transmitter frequencies 
of 600.1517482 MHz (^H) and 242.9458642 MHz (^^P). Ti re- 
laxation times were measured by inversion-recovery methods, 
while T2 times were measured using a single spin-echo; all 
relaxation times are averaged over the four components of 
the relevant multiplet, as there was little sign of correlated 
relaxation. 
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Na2C03 to remove the silyl protecting group ||22l. HCl 
was added cis across the triple bond using LiCl in acetic 
acid p7t , and the Z-E isonierisation was accomplished by 
heating with PhSH in the presence of AIBN |^ . Acidic 
hydrolysis of the ester groups gave E-(2-chloroethenyl)- 
phosphonic acid (overall yield 5%, purity > 95%), from 
which contaminating metal ions were removed using 
chelex resin. The NMR sample was prepared by dis- 
solving 7mg of compound into 600/iL of D2O to give a 
0.08M solution, and all NMR experiments were run at a 
nominal temperature of 20° using a homebuilt 600 MHz 
(^H frequency) NMR spectrometer at the OCMS with a 
homebuilt double resonance (^H inner, ^^P outer) probe. 
The measured NMR parameters are listed in Table ^ 

This system may be conveniently described using prod- 
uct operator notation p^ . The experiment begins with 
the preparation of the initial state PzAqBq using a mod- 
ification of the "cat state" method of Knill et al. [ p2[ . 
Qubit P can then be set to any desired point on the 
Bloch sphere using a single radiofrequency (RF) pulse. 
After the cloning sequence, the reduced density matrices 
of qubits A and B will correspond to Bloch vectors ly- 
ing parallel to the original Bloch vector of P, but with 
lengths only 2/3 that of the original vector; note that the 
effect of the maximally mixed component is simply to re- 
duce the length of the Bloch vector without affecting its 
orientation. Detailed pulse sequences are shown in Fig. ^ 
and described below. 

NMR pulse sequences were developed by replacing the 
abstract gates in Fig. |^ with an idealized sequence of 
NMR pulses (including 2-rotations) and delays; the re- 
sulting sequences were then simplified by adsorbing z- 
rotations into abstract reference frames |2^, Q and com- 
bining RF pulses when convenient |2^] . The detailed im- 
plementation of the sequence was chosen so as to avoid 
selective pulses as far as possible; those selective pulses 
in the "echo" sequence that could not be avoided were 
implemented using a variant on the "jump and return" 
sequence 18, This sequence incorporates a delay 

which should naively take the value ego = \/{A5v), where 
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FIG. 3: NMR pulse sequences. White and black boxes are 
90° and 180° pulses, while grey boxes are pulses with other 
flip angles; pulse phases and gradient directions are shown be- 
low each pulse. The large ovals correspond to the "echo" se- 
quence described below. All RF pulses are hard, with ^H fre- 

fr icy selection achieved using "jump and return" sequences 
M. The purification sequence (a) and cloning sequence 
are built around the echo sequence (c) which implements 
the coupling element of PA and PB controlled phase shifts 
pif . Delay times are tab = 1/(4: Jab), tap = 1/(8 Jap) and 
tbp — l/(8Jsp); ego was set as described in the main text. 



the two signals have frequencies -ii5v, but in practice 
it is better to slightly reduce its length in order to al- 
low for non- idealities arising from the fact that 5v is not 
much greater than the J couplings; in our experiments 
ego was reduced to 90% of its nominal value. 

The initial "purification" sequence was developed sep- 
arately, but uses many of the same ideas. Our implemen- 
tation is based on the "cat state" methods |2^ of Knill 
et al. The core of the sequence comprises a Hadamard 
gate and a pair of controUed-NOT gates that convert 
to a density matrix containing all terms along the anti- 
diagonal, a gradient filter that selects the desired triple 
quantum term, and a second pair of controUed-NOT gates 
and a Hadamard which act to produce PzAqBq. The 
novel gradient filter sequence used here uses two gradi- 
ents with strengths in the ratio 1 : 0.6633 to select the 
desired three quantum terms while crushing other anti- 
diagonal terms; note that phase cycling is not required. 

Unlike conventional gradient sequences used to select 
coherence transfer pathways ||, ^ , this scheme is largely 
unaffected by diffusion. It is, however, vulnerable to im- 
perfections in the anticat sequence: these are removed by 
the final z-filter, which incorporates a variable delay to 
suppress zero quantum terms |T^. Note that this vari- 
able delay can be combined with the standard CYCLOPS 
phase cycle Q , and so no additional phase cycling is nec- 
essary ]l8| . In principle this sequence will also suppress 
terms arising from initial Az and B^ magnetization, as 
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FIG. 4: The experimental result from cloning the initial state 
Px. The receiver phase was set using a separate experiment 
so that x-magnetization appears as positive absorption mode 
lines. The NMR experiment was repeated four times using the 
CYCLOPS phase cycle H to reduce instrumental imperfections. 



well as signal from residual HOD and other impurities; 
in reality the behaviour is not as selective as one might 
desire, and so the purification sequence is preceded by 
a 90° excitation pulse and gradient crush, to reduce 
any signals that do not start from P^. 

The effects of the idealized network (Fig. |]) on an NMR 
system with spin P set to an arbitrary point on the Bloch 
sphere, Pg^, are easily calculated: the observable NMR 
signal takes the form 
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[Aec^ {BaPo + BiPi) + Be^ {A^Pa + ^iPi)] . (1) 



After tracing out the other two spins the signal from 
spins A and B is reduced to the expected forms, 
and '^Bg^ respectively. Note that the component corre- 
sponding to the maximally mixed state has no effect on 
the NMR signal as all NMR observables are traceless. 
Spin P could be traced out experimentally by applying 
a '^^P decoupling field during observation of the two 
nuclei, A and B, but this approach was found to lead 
to unacceptable sample heating, with resultant shifts in 
resonance frequencies. Similarly it would in principle be 
possible to trace out spin A while observing B by the use 
of selective homonuclear decoupling, but this has the fur- 
ther disadvantage of requiring two separate experiments 
to observe spins A and B. For these reasons we decided to 
implement the tracing out process in software by adding 
together signals arising from the four components in each 
multiplet; this is most simply achieved by integrating the 
entire multiplet. 

The experimental results from cloning the initial state 
Px are shown in Fig. ^. The spectrum does have the 
expected overall form (inphase absorption signals at the 
outermost positions of each multiplet), but some non- 
idealities are clearly visible: the two lines in each multi- 
plet have different heights, and there are small signals at 
other positions in the multiplets. More seriously the ob- 
served signal intensity is only about 25% of the expected 
value. Detailed calculations (data not shown) indicate 
that the errors in relative line intensities can be traced 
to J coupling evolution during selective excitation se- 
quences, while the overall intensity loss can be ascribed 
to T2 relaxation and Bi inhomogeneity effects. 



An important feature of the approximate cloning net- 
work is that all input states are cloned equally well, and 
so it is necessary to study the behaviour of the pulse 
sequence when applied to a wide range of points on the 
Bloch sphere. We have studied a total of 312 input states, 
arranged on a 13 by 24 rectangular grid of 9 and val- 
ues, with a spacing of 15°: this choice does not cover the 
Bloch sphere uniformly (for example the north and south 
poles are both sampled 24 times), but is experimentally 
convenient. For each input state we measured the total 
NMR signals observed from spins A and B, and their real 
and imaginary components are plotted in Fig. |^. The ex- 
perimental results clearly show the expected cosine and 
sine modulations, indicating that the cloning network is 
effective for all these input states. 

Finally, we recall that approximate quantum cloning 
is unitarily equivalent to another interesting operation, 
the optimal universal-NOT gate and so our exper- 
iment could also be considered as an implementation 
of universal-NOT; this point is not, however, specifically 
demonstrated here. 
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integrals of the real (a, b, c) and imaginary (d, e, f) parts of the NMR signals from spin A (a, d) and spin B (c, f) are shown 
as mesh and contour plots as a function of 9 and 0; the theoretical results are shown in (b, e). 
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